Shear-current effect in a turbulent convection with a large-scale shear 
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The shear-current effect in a nonrotating homogeneous turbulent convection with a large-scale 
constant shear is studied. The large-scale velocity shear causes anisotropy of turbulent convection, 
which produces the mean electromotive force £^- w ^ oc WxJ and the mean electric current along 
the original mean magnetic field, where W is the background mean vorticity due to the shear and 
J is the mean electric current. This results in a large-scale dynamo even in a nonrotating and 
nonhelical homogeneous sheared turbulent convection, whereby the a effect vanishes. It is found 
that turbulent convection promotes the shear-current dynamo instability, i.e., the heat flux causes 
positive contribution to the shear-current effect. However, there is no dynamo action due to the 
shear-current effect for small hydrodynamic and magnetic Reynolds numbers even in a turbulent 
convection, if the spatial scaling for the turbulent correlation time is r(k) oc k~ 2 , where k is the 
small-scale wave number. 

PACS numbers: 47.65.Md; 47.65.-d 



INTRODUCTION 



It is generally believed that the large-scale magnetic 
fields of the Sun, solar type stars and galaxies are origi- 
nated by a dynamo process, i.e., due to a joint action of 
small-scale helical turbulent motions (the a effect) and 
large-scale differential rotation (see [p. B, j, B 1 I \ 7 , B, 
9]). It has been recently recognized flQ. Illl.ll3ll3l.ll4 Il9] 
that in a sheared nonhelical and nonrotating homoge- 
neous turbulence whereby the a effect vanishes, the 
mean-field dynamo is possible due to shear-current ef- 
fect. 

The mechanism of the shear-current dynamo is as fol- 
lows. The deformations of the original nonuniform mag- 
netic field lines are caused by upward and downward tur- 
bulent eddies. In a sheared turbulence the inhomogeneity 
of the original mean magnetic field breaks a symmetry 
between the influence of the upward and downward tur- 
bulent eddies on the mean magnetic field. This creates 
the mean electric current J along the mean magnetic field 
B and produces the shear-current dynamo. In particu- 
lar, the large-scale velocity shear creates anisotropy of 
turbulence. This produces the mean electromotive force 
S^ w%> oc WxJ, where W is the background mean vortic- 
ity due to the shear. Joint effects of the mean electromo- 
tive force £( w ) and stretching of the mean magnetic field 
due to the large-scale shear motions cause the mean-field 
dynamo instability. 

A sheared turbulence is a universal feature in astro- 
physics. The shear-current effect might be an origin for 
the large-scale magnetic fields in colliding protogalactic 
clouds and in merging protostellar clouds [l5|. This ef- 
fect might be also important in accretion discs where 
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the mean velocity shear comes together with rotation, so 
that both the shear-current effect and the a effect might 
operate. Since the shear-current effect is not quenched 
(see [HI, [3]) contrary to the quenching of the nonlinear 
a effect, the shear-current effect might be the only sur- 
viving effect, and it can explain the origin of large-scale 
magnetic fields in astrophysical plasmas with large-scale 
sheared motions. 



The shear-current effect is a fundamental phenomenon 
which should be studied in different situations, e.g., in 
different types of turbulence. The goal of the present 
study is to investigate the shear-current effect in a non- 
rotating homogeneous turbulent convection with a large- 
scale constant velocity shear. Note also that in many 
astrophysical applications turbulent convection plays an 
important role, e.g., in the convective zones of the Sun 
and solar type stars. We have shown that turbulent con- 
vection promotes the shear-current dynamo instability. 
In particular, the heat flux causes positive contribution 
to the shear-current effect. However, the shear-current 
dynamo is impossible for small hydrodynamic or mag- 
netic Reynolds numbers even in a turbulent convection, 
if the spatial scaling for the turbulent correlation time is 
r(k) oc k~ 2 , where k is the small-scale wave number. 



This paper is organized as follows. In section II we 
formulate the governing equations, the assumptions and 
the procedure of the derivations. In section III we study 
properties of the shear-current effect in a sheared tur- 
bulent convection and discuss the shear-current dynamo. 
In section IV we draw concluding remarks and perform 
comparison of the theoretical predictions with the direct 
numerical simulations. Finally, in Appendix A we per- 
form a detailed derivation of the shear-current effect in a 
turbulent convection. 
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II. THE GOVERNING EQUATIONS 

In order to study the shear-current effect in a turbu- 
lent convection we use a procedure which is similar to 
that applied in [IHtnl- I n particular, we employ a mean 
field approach whereby the pressure, entropy, velocity 
and magnetic fields are separated into the mean and fluc- 
tuating parts, where the fluctuating parts have zero mean 
values. To determine the effect of shear on a turbulent 
convection we use equations for fluctuations of velocity, 
magnetic field and entropy 



^ = _(U-V)u-(u-V)U-v(£ 

at \p 

+ — [(b-V)B + (B-V)b]+u JV , 

Po 



gs 



db 



— = (B-V)u-(u-V)B + (b-V)U-(U-V)b + b JV , 



ds 
dt 
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where B, U and S are the mean magnetic field, the mean 
velocity field and the mean entropy, u, b and s are fluc- 
tuations of velocity, magnetic field and entropy, po is the 
fluid density, TV 2 = — g- and g is the acceleration of 
gravity, e is the unit vector directed opposite to g, the 
magnetic permeability of the fluid is included in the def- 
inition of the magnetic field, p are the fluctuations of to- 
tal (hydrodynamic and magnetic) pressure, v", h N and 
s N are the nonlinear terms which include the molecular 
viscous and diffusion terms. Equations U)-© for fluc- 
tuations of fluid velocity, entropy and magnetic field are 
written in the Boussinesq approximation. We consider 
the hydrostatic nearly isentropic basic reference state. 
The turbulent convection is regarded as a small devia- 
tion from a well-mixed adiabatic reference state. 

Using Eqs. H])-© written in a Fourier space we derive 
equations for the instantaneous two-point second-order 
correlation functions of the velocity fluctuations (utUj), 
the magnetic fluctuations (bi bj), the entropy fluctuations 
(ss), the cross- helicity tensor (biUj), the turbulent heat 
flux {sui) and (sbi). The equations for these correlation 
functions are given by Eqs. (|A.2[) - (|A.7p in Appendix A. 
We split the tensor of magnetic fluctuations into non- 

( H ) 

helical, hij = (bibj), and helical, hh , parts. The helical 

part h\j depends on the magnetic helicity, and it is de- 
termined by the dynamic equation which follows from 
the ma gne tic helicity conservation arguments (see, e.g., 
[3 HM SI II M H H! , and a review @ ) . 

The second-moment equations include the first-order 
spatial differential operators N applied to the third-order 
moments M {III \ A problem arises how to close the 
system, i.e., how to express the set of the third-order 
terms SfM^- 111 ^ through the lower moments M^- 11 ^ (see, 
e.g., [H, [13, We use the spectral r approxima- 

tion which postulates that the deviations of the third- 



moment terms, AfM^ 111 " 1 (k), from the contributions to 
these terms afforded by the background turbulent con- 
vection, AfM^ 111 ' ^ (k), are expressed through the similar 
deviations of the second moments, M^ // )(k) — M( 77 '°)(k): 

I 



r(fc) 



[M( /J) (k)- M (// ' 0) (k)] , (4) 



(see [1J, M, M, [20, Hll), where r(fe) is the scale- 
dependent relaxation time, which can be identified with 
the correlation time of the turbulent velocity field. The 
quantities with the superscript (0) correspond to the 
background shear-free turbulent convection with a zero 
mean magnetic field. We apply the spectral r approxi- 
(1) mation only for the nonhelical part of the tensor of 
magnetic fluctuations. Note that a justification of the 
t approximation for different situations has been per- 
formed in numerical simulations and analytical studies 

in d, us, m, m he hi . 

We assume that the characteristic time of variation 
of the mean magnetic field B is substantially larger 
than the correlation time r(fc) for all turbulence scales. 
This allows us to get a stationary solution for the equa- 
tions for the second-order moments, M^ lr \ We split 
all second-order correlation functions, M^ 11 ', into sym- 



metric h 



h 



(a) 



ij 



[Tiy(k) + hij(— k)]/2 and antisymmetric 

[/ijj(k) — hij(— k)]/2 parts with respect to the 
wave vector k. For the integration in k-space we have to 
specify a model for the background shear-free turbulent 
convection (i.e., a turbulent convection with B = 0). The 
background turbulent convection is maintained by an im- 
posed vertical heat flux F* = (su z ) with divF* = at 
a low boundary of convective region. We used the fol- 
lowing model for the homogeneous background turbulent 
convection: 



K^) (0) (k) 
(W (0) (k) 

(k) 



■:u-)Ptf(*)W(fc) 



2 > 

b 2 )P y -(fc)IU(£0 , 



,(o) 



3 (s u z ) e r , 
kj I k ^ 



,(*)w(*o 



(5) 
(6) 
(7) 



6ij is the Kronecker 



where Pij(k) = 5, 

tensor, W(k) = E(k)/8irk 2 , the energy spectrum is 
E(k) = (q- l)(k/k )-?, k = l/l and the length l 
is the maximum scale of turbulent motions. The tur- 
bulent correlation time is r(fc) = Cto (/c//co) _M , where 
the coefficient C — (q — 1 + p.)/(q — I). This value of 
the coefficient C corresponds to the standard form of the 
turbulent diffusion coefficient in the isotropic case, i.e., 
Vt = I T ( k ) [(u 2 ) E(k)]dk = t (u 2 )/3. Here the time 
To = fo/V (u 2 ) and \J (u 2 ) is the characteristic turbu- 
lent velocity in the scale 1$. For the Kolmogorov's type 
background turbulence (i.e., for a turbulence with a con- 
stant energy flux over the spectrum), the energy spec- 
trum E(k) oc —dr/dk, the exponent fx = q — 1 and the 
coefficient C = 2. In the case of a turbulence with a scale- 
independent correlation time, the exponent p, = and 
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the coefficient C = 1. Motions in the background tur- 
bulent convection are assumed to be non-helical. Using 
the solution of the derived second- moment equations, we 
determine the contributions to the mean electromotive 
force, £[ — Ei mn J (b n «m)k dk, caused by the sheared 
turbulence (see Appendix A), where Sijk is the fully an- 
tisymmetric Levi-Civita tensor. This procedure allows us 
to determine the contributions to the shear-current effect 
caused by the sheared turbulent convection. 



III. THE SHEAR-CURRENT DYNAMO 

We consider a homogeneous turbulent convection with 
a constant mean velocity shear, U = (0, 5a;, 0) and W = 
(0,0, S) . We consider a most simple form of the mean 
magnetic field, B = (B x (z), B y (z), 0). The contributions 
to the mean electromotive force caused by the sheared 



by 



turbulence, are £° 



bf jk is given by 



bf jk Vfc Bj, where the tensor bf jk 



MVZ7) 2 + ry r (Vi?) 2 /p 

g{su z ) 



(10) 



where v T is the turbulent viscosity and r] T is the coeffi- 
cient of turbulent magnetic diffusion. 

Therefore, in the kinematic approximation the mean 
magnetic field is determined by 



dB x 
dt 

dBy 

dt 



-a B Sl 2 B" + ri T BZ, 



S B x 



l 



(11) 

(12) 



where B" = d 2 Bi/dz 2 . Here we neglect small contri- 
butions to the coefficient of turbulent magnetic diffusion 
caused by the shear motions because Stq <C 1. The di- 
mensionless parameter a B describes the shear-current ef- 
fect. Straightforward calculations using Eqs. (JH])-© yield 
the parameter a B — a™ + er^, where 



°ijk 



l 2 o ^ dkn [Qo V n Uj + 2 Q 1 [dU) nj ] , (8) 

,2 h 

a * l ° 140 



Uj + [Q 2 e lkn 



~\~Q4 i^ikm &mn ^inrn ^mfc)] (^^0 



"J 



(9) 



Equations ([5]) and © are derived in Appendix A. Here 
(dU)ij = {ViUj + V ' jUi)/2 and the coefficients Q„ are 
Qo = (3-2 fi)-e (5+2 /x), Qi = e (7+6 Q 2 = 

Q3 = 18 — 19/j,, Qi = /1 — 6, the parameter e = E m / E v , 
E m and E v are the magnetic and kinetic energies per 
unit mass in the background turbulent convection, and 



h = / T 2 (k)E(k)dk 

h / T 3 (k)E(k)dk 



(9-1 + M) 2 



{q - l + 2 M ) (q - 1) 

(g - 1 + M) 3 
(q - 1 + 3 M ) (q - l) 2 



For the Kolmogorov's type turbulence, the exponent 
fi = q — 1 and the parameters 7 2 = 4/3 and I3 = 2. 
In the case of a turbulence with a scale-independent cor- 
relation time, the exponent /i — and the parameters 
1-2 = I3 = 1. The tensor b^ k in Eq. © describes the 
contributions of the heat flux to the shear-current effect, 
while tensor b^ k determines the non-convective contri- 
butions (which are independent of the heat flux) to the 
shear-current effect. In Eqs. ([S])-® we have taken into 
account only the terms which contribute to the shear- 
current effect. In particular, we have taken into account 
that By B x (see Eq. (|2"0|) below) and considered a 
weak mean velocity shear U = (0, Sx, 0), where St <C 1. 
The convective contribution to the dynamo instability 
due to the shear-current effect depends on the parame- 
ter a» = 2cjtqF* I (u 2 ) which is determined by the budget 
equation for the total energy. The parameter a* is given 



30 



(Qo + Qi 



280 



Q 2 - Qi + 2 (Q 3 + Qi) sin 2 



(13) 



(14) 



4> is the angle between the unit vector e and the back- 
ground vorticity W due to the large-scale shear. Equa- 
tions (fT3|) and (fT4|) yield the following final expressions 
for the parameter a B : 



h 
15 



1- //-:- mi 2 M ) + ^ 3 -[2 



7/2 



+3 (2 - 3/x) sin 2 cf>] 



(15) 



where the terms cx a* in Eq. (|15p describe the contri- 
bution of the turbulent convection to the shear-current 
effect. Equations (fTTj) and (TT2"|) determine the shear- 
current dynamo instability. In particular, the first term 
cx SB X in the right hand side of Eq. (JT2J) determines 
the stretching of the magnetic field B x by the shear mo- 
tions and produces the field B y . On the other hand, the 
interaction of the non-uniform magnetic field B y with 
the background vorticity W produces the electric cur- 
rent along the field B y . This effect is determined by the 
first term in the right hand side of Eq. (jTT|) and causes 
the generation of the magnetic field component B x . The 
growth rate of the mean magnetic field due to the shear- 
current dynamo instability is given by 



y = Sl y/o^Kz - Vt K z 



(16) 



where K z is the large-scale wave number. The necessary 
condition for the dynamo instability is a B > 0. 

The shear-current dynamo instability depends on the 
spatial scaling of the correlation time r(k) cx k~^ of 
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the turbulent velocity field, where k is the small-scale 
wave number. In the absence of turbulent convection, 
the terms oc a* in Eq. (|15p vanish, and the shear-current 
dynamo in a non-convective turbulence with e = oc- 
curs for fj, < 1. For the Kolmogorov's type turbulence, 
the exponent fi — 2/3 and Eq. (fT5|) reads 



4 

135 



1 + 7e 



(17) 



In this case the parameter a B is independent of the angle 
<f> between the unit vector e and the background vorticity 
W. For a turbulent convection with a scale-independent 
correlation time, the exponent \x = and the parameter 
a B is given by 



1 

15 



1 



e + — a* (1 + 3 sin 



(18) 



In these cases the shear-current dynamo instability 
causes the generation of the large-scale magnetic field. 
It follows from Eqs. (fl"5|) - (fT5|) that turbulent convection 
promotes the shear-current dynamo instability. In par- 
ticular, the heat flux causes positive contribution to the 
shear-current effect when 2 + 3 (2 — 3/i) sin 2 > 0. 

However, for small hydrodynamic and magnetic 
Reynolds numbers, the turbulent correlation time is of 
the order of r(fc) oc l/(vk 2 ) or r(fc) oc l/(r/fc 2 ) depending 
on the magnetic Prandtl number, i.e., r(k) oc fc~ 2 . In 
this case fi = 2, and the parameter a B < even in a tur- 
bulent convection with e = 0. This implies that for small 
hydrodynamic and magnetic Reynolds numbers there is 
no dynamo action due to the shear-current effect. This 
result is in agreement with the recent studies [37ll3l| per- 
formed in the framework of the second order correlation 
approximation (SOCA) for sheared non-convective flows. 
This approximation is valid only for small hydrodynamic 
Reynolds numbers. Even in a highly conductivity limit 
(large magnetic Reynolds numbers), SOCA can be valid 
only for small Strouhal numbers, while for large hydro- 
dynamic Reynolds numbers (for a developed turbulence), 
the Strouhal number is 1. 

Note that the standard approach (i.e., SOCA) cannot 
describe the situation in principle. The reason is that the 
shear-current dynamo requires a finite correlation time 
of turbulent velocity field, so the delta-correlated version 
of SOCA fails. The application of the path integral ap- 
proach for the study of the shear-current dynamo also re- 
quires a finite correlation time of turbulent velocity field. 
The shear-current dynamo is a phenomenon that results 
from the interaction of the energy-containing-scale of 
turbulence with large-scale shear, and the constraint is 
that the hydrodynamic and magnetic Reynolds numbers 
should be not small at least. Therefore, the SOCA-based 
approaches do not work properly to describe the shear- 
current dynamo. Probably, the hydrodynamic and mag- 
netic Reynolds numbers can be of the order of unity and 
there is no need for a developed inertial range in order 
to maintain the shear-current dynamo. 



In order to determine the threshold required for the ex- 
citation of the shear-current dynamo instability, we con- 
sider the solution of Eqs. (jTTJ) and (TT2")) with the following 
boundary conditions: B(t, \z\ = L) = for a layer of the 
thickness 2L in the z direction. The solution for the mean 
magnetic field is determined by 



B y (t,z) = Bq exp(7*) cos(K z z + p) , 



B x {t,z) = l K z 



■B y (t,z) 



(19) 
(20) 



For the symmetric mode the angle p = tt n and the 
large-scale wave number K z — (7t/2)(2to + 1) L -1 , where 
n, m = 0, 1, 2, ... . For this mode the mean magnetic field 
is symmetric relative to the middle plane z = 0. Let us 
introduce the dynamo number D = (Iq S*/L) 2 a B , where 
parameter S* — S L 2 /rj T is the dimensionless shear num- 
ber. For the symmetric mode the mean magnetic field 
is generated due to the shear-current effect when the dy- 
namo number D > D CI = {it 2 / '4) (2m + l) 2 . For the 
antisymmetric mode the angle tp — (ir/2) (2n + 1) with 
n = 0, 1, 2, the wave number K z = nmL^ 1 and the 
magnetic field is generated when the dynamo number 
D > D CI = TT 2 m 2 , where m = 1,2,3,.... The max- 
imum growth rate of the mean magnetic field in the 



shear-current dynamo instability, j n 



S 2 l 2 



is attained at K z = SIq ^Ja B j2r\ T . Therefore, the char- 
acteristic scale of the mean magnetic field variations 
Lb = 2n/K z — &uq/(S ^fa^) . For the shear-current dy- 
namo, the ratio of the field components B x /B y is small 
(see Eq. (f20| ). Remarkably, in the afl dynamo, the 
poloidal component of the mean magnetic field is much 
smaller than the toroidal field. 



IV. DISCUSSION 

In the present study we investigate the shear-current 
effect in a nonrotating homogeneous turbulent convec- 
tion with a large-scale constant velocity shear. We show 
that the condition for the shear-current dynamo is in- 
dependent of the exponent of the energy spectrum of 
turbulent convection, but it depends on the scaling ex- 
ponent fi of the turbulent correlation time r(k) oc fc~ M , 
where k is the small-scale wave number. We discuss three 
cases in details: (i) the Kolmogorov's type turbulence 
with the exponent fi = 2/3; (ii) a turbulent convection 
with a scale-independent correlation time (/z = 0); (iii) a 
turbulent convection with small hydrodynamic and mag- 
netic Reynolds numbers with the scaling r(fc) oc k~ 2 . We 
have found that turbulent convection promotes the shear- 
current dynamo instability. In particular, the heat flux 
causes positive contribution to the shear-current instabil- 
ity. However, the shear-current dynamo does not occur 
for small hydrodynamic and magnetic Reynolds numbers 
even in a turbulent convection, if the spatial scaling for 
the turbulent correlation time is r(k) oc k~ 2 . 

For simplicity we consider weak linear velocity shear, 
U = (0,Sx,0), where the parameter St <C 1. The 
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main effect of the weak linear velocity shear on turbu- 
lent convection is a generation of additional anisotropic 
velocity fluctuations. We consider turbulent convection 
in the region which is far from the boundaries, because 
the constant linear velocity shear cannot exist near the 
boundaries whereby the boundary layers form. The gen- 
eration of the magnetic field in a nonlinear velocity shear 
depends on boundary conditions and requires numerical 
study. Turbulent convection can be inhomogeneous in 
this case. 

The main goal of this paper is to study an effect of 
the heat flux on the shear-current dynamo instability in 
a most simple model of turbulent convection with a lin- 
ear shear. The shear-current dynamo acts also in inho- 
mogeneous turbulent convection. However, in inhomo- 
geneous turbulence with a large-scale constant velocity 
shear the kinetic helicity and the a effect do not vanish 
(see [13, HUES). 

In this case the shear-current dynamo 
acts together with the a-shear dynamo which is similar 
to the afl dynamo. The joint action of the shear-current 
and the a-shear dynamos have been recently discussed in 
[l6L I39L l4 0l|. The shear-current effect does not quenched 
(see [111 . Il6| ] ) contrary to the quenching of the nonlinear a 
effect, the turbulent magnetic diffusion, the effective drift 
velocity, etc. Therefore, the shear-current effect might be 
the only surviving effect, and it can explain the origin of 
large-scale magnetic fields in sheared astrophysical tur- 
bulence. 

The shear-current dynamo instability is saturated by 
the nonlinear effects. The nonlinear mean-field dynamo 
due to a shear-current effect in a nonhelical homoge- 
neous turbulence with a mean velocity shear has been 
investigated recently in [3] (see also [40l|). whereby the 
transport of magnetic helicity as a dynamical nonlinear- 
ity has been taken into account. The magnetic helicity 
flux strongly affects the saturated level of the mean mag- 
netic field in the nonlinear stage of the dynamo action. 
In particular, numerical solutions fl4j ] of the nonlinear 
mean-field dynamo equations which take into account the 
shear-current effect, show that if the divergence of the 
magnetic helicity flux is not small, the saturated level of 
the mean magnetic field is of the order of the equipar- 
tition field determined by the turbulent kinetic energy. 
These results are in a good agreement with direct nu- 
merical simulations [H, [HI], whereby the generation of 
the large-scale magnetic field in a nonhelical turbulence 
with an imposed mean velocity shear has been investi- 
gated. 

In the direct numerical simulations El IT3I ] the non- 
convective turbulence is driven by a forcing that consists 
of eigenfunctions of the curl operator with the wavenum- 
bers 4.5 < kf < 5.5 and of large-scale component with 
wavenumber k\ = 1 . The forcing produces the mean flow 
U = Uo cos (fci x) cos (ki x). The numerical resolution in 
these simulations is 128 x 512 x 128 meshpoints, and the 
parameters used in these simulations are as following: 
the magnetic Reynolds number Rm = U Tms /(r] kf) = 80, 
the magnetic Prandtl number Pr m = v/rj = 1 and 



Uo/u Tms — 5. The growth rate of the mean magnetic 
field is about 7 To w 2 x 1CP 2 . This allows us to estimate 
the parameter a B characterizing the shear-current effect, 
a B rs 3.3 x 10~ 2 . On the other hand, our theory predicts 
(i B = (3 — 6) x 10~ 2 depending on the parameter fi. Note 
that in the numerical simulations [H, EH the- shear is 
not small (i.e., the parameter Stq ~ 1), which explains 
some difference between the theoretical predictions and 
numerical simulations. Therefore, the numerical simu- 
lations [H EH clearly demonstrate the existence of the 
large-scale dynamo in the absence of mean kinetic helic- 
ity and alpha effect, in agreement with the theoretical 
predictions discussed in the present paper. 
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APPENDIX A: THE ELECTROMOTIVE FORCE 
IN A SHEARED TURBULENT CONVECTION 

In order to study the shear-current effect in a sheared 
turbulent convection we use a procedure applied in (Til . 
fl7j for similar problems. Let us derive equations for the 
second moments. To exclude the pressure term from the 
equation of motion ([I]) we calculate Vx(Vxu). Then 
we rewrite the obtained equation and Eqs. ©-([H]) in a 
Fourier space. We also apply the two-scale approach, 
e.g., we use large scale R = (x + y)/2, K = ki + k 2 and 
small scale r = x — y, k = (ki — k 2 )/2 variables (see, 
e.g., [H|). This implies that we assume that there exists a 
separation of scales, i.e., the maximum scale of turbulent 
motions Iq is much smaller then the characteristic scale 
Lb of inhomogeneities of the mean magnetic field. We 
derive equations for the following correlation functions: 

/ij(k) = L(ui\Uj) , hij(k) = L(bi-,bj) , 

3tj(k) = L(bi;uj) , Fj(k) = L(s; Ui) , 

Gi(k) = L(s;h) , 9(k) = L(s;s) , (A.l) 

where 

L(a; c) = J (o(k + K/2)c(-k + K/2)) 

x exp (iK-R) dK . 

The equations for these correlation functions are given 
by 

^jP- = i(k-B)$« + 4 + I? jmn (U)f mn +Aff ij , 
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(A.2) 



dhjj (k) 
dt 

%j( k ) 
dt 



a* 

ag<(k) 

dt 

OQ (k) 
dt 



(k-B)^ + + ^ mn (U)/i m „ + A/7ij 



(A.3) 



i(k-B)[/«(k)-/ Hj -(k)-/iW] + /g 
+^m™( U )-9™« + ffe n P J „(fc)G i (-k) 

+^<% , (A.4) 
-i(k.B)G«(k)+ if + /ff m (U)F m 

+fle„P in (fe)e(k)+^F i , (A.5) 
-i (k-B)Fi(k) + if + (V m l/,) G m (k) 

(A.6) 

-fiCkj+^e, (A.7) 



N 2 
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where hereafter we omit argument t and R in the corre- 
lation functions and neglect small terms ~ 0(V 2 ). Here 
V = 8/ OR, and we also neglect a small term oc N 2 /g in 
Eq. (33). InEqs. $„(k) = fe (k)- 5ji (-k), 

Py(fc) = % - kikjjk 2 , Affij_ = ge n [P m {k)F {k) + 
Pj n (k)Fi(-k)] + Jiffy, and Aff ij; Mh tJ , tf gij , AfF t , 
NGi and JvQ are the third-order moment terms appear- 
ing due to the nonlinear terms. The terms which are 
proportional to the heat flux F{ in the tensor Affij , can 
be considered as a stirring force for the turbulent convec- 
tion. Note that a stirring force in the Navier-Stokes tur- 
bulence is an external parameter. The tensors IZ mn (U), 
Etjmn (U) , J^ m „ (U) and H% (U) are given by 



I a (U) 

7 ruin V ^ / 



1kiq8 7n p8 j n -f" QikjqSimSpn ^im^jq^pn 

d 



t (U) 



$im$jq3pn ~\~ ^jm^iq^pn 

d 



dk r . 



V P U q , 



~\~Si m Sj n kq 



tjmn \ J 



v P c/ 9 , 



where fcjj = kikj/k 2 . The source terms , i|j, iy, 
if and If which contain the large-scale spatial deriva- 
tives of the mean magnetic field, are given in [ill \v\ 
(see also Eqs. (|A.12j) and (|A.13|) below). Next, in 
Eqs. (|A.2p - (|A.7[) we split the tensor for magnetic flue- 
tuations into nonhelical, hjj, and helical, , parts. 
The helical part of the tensor of magnetic fluctuations 



h\j depends on the magnetic helicity and it follows from 
the magnetic helicity conservation arguments (see, e.g., 
[3, 0, S3, SI S3, SI)- We also use the spectral r 
approximation which postulates that the deviations of 
the third-moment terms, NM^ 111 ' (k) , from the contri- 
butions to these terms afforded by the background tur- 
bulent convection, UM^ nl ^(k), are expressed through 
the similar deviations of the second moments, M^ 11 * 1 (k) — 
Af( /7 <°)(k) [see Eq. ©]. 

We take into account that the characteristic time 
of variation of the mean magnetic field B is substan- 
tially larger than the correlation time r(k) for all tur- 
bulence scales. This allows us to get a stationary solu- 
tion for Eqs. (|A.2[) - (|A.7[) f° r the second-order moments, 
A/( /7 )(k), which are the sum of contributions caused by a 
shear-free turbulent convection and a sheared turbulent 
convection. The contributions to the mean electromo- 
tive force caused by a shear-free turbulent convection are 
given in [l7| ■ On the other hand, the contributions to the 
mean electromotive force caused by the sheared turbu- 
lent convection are £^ = £ m ji J <7.y(k) dk. In particular, 
in the kinematic approximation the contributions to the 
cross-helicity tensor g^ caused by the sheared turbulent 
convection, are given by 



9ijW — T [Jijmn9mn 



■ge n P jn (k)G?(-K)] , 
(A. 



where 



Gf(k) 

9ij 



= T 2 (V m U t )lf , (A.9) 
= T[I? j +rge n P jn (k)lf] , (A. 10) 

= T[(2P js (k)-6 js )E? kmn h<$ n 



A T a f(0)| d 
°»s 1 kjmn Jmnl ^Sjk 

&ii $mk ~t~ 77 &im kj 



dk k 



(All) 
(A.12) 



if. 

i j 



(2P jn (k)-S jn )h^-8 in fif 



2 " dk k U « + « ' 



B n , k , (A.13) 



and Bij — V jBi. We take into account that in Eq. (|A.8[) 
the terms with symmetric tensors with respect to the in- 
dexes "i" and "j" do not contribute to the mean elec- 
tromotive force because £ ^ = e m ji J gfj (k) dk. For 
the integration in k-space we have to specify a model 
for the background shear-free turbulent convection (with 
B = 0), which is determined by Eqs. I©-© in Section 
III. 

The contributions to the mean electromotive force 
caused by the sheared turbulent convection, are £ f = 



b ijk V fc B j> where the tensor b ijk = b ijk ^ 
Eqs. ©-© in Section III. For derivation of Eqs. ©-© 
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we use the following identities: 



k^k^k ra k ri . 
— sin dO dip 



k 4 



k* ■ L? ■ Ic L? k* Ic 

- sin 6 dO dip 



k 6 



and 



4tt 
15 
4tt 



A, 



105 



A,- 



In Eqs. ([8])- ([9]) we have taken into account only the terms 
which contribute to the shear-current effect. In particu- 
lar, we consider the mean magnetic field in a most simple 
form B = (B x (z), B y (z), 0) and we take into account that 
By 3> B x and S tq -C 1, where the mean velocity shear 
is U = (0, Sx, 0) and W = (0, 0, S). Straightforward cal- 
culations using Eqs. ©-© and Eqs. (13]) and Q3|) yield 
<r„ = <j u „ + , where 



^ijmnpq ^mnpq &ij 



^jmnq &ip 



^imnq &jp 



-A, 



>Si q 



and 



^ p Uq — 2 £ikr 



[(dU) 



mj 



2 €qj Vm Uq\ , 



H~2 Gmp L/j H" 
£inm &ms ^jkpqns ^p Uq — 4 {iSmm [Cmfc \^Ujjij 

(dU) kn }} 

£inq Cms ^jkpmns ^p Uq — &inq \^ p Uq [2 Gr^ft /S.jp mn 

+2e r 



~-mj ^kpmn ~^ A^pfc n ] V n £/g (^jfc "I - ^Gjfc)} , 



£?jn ^ttjs ^kpqmns ^p Uq — 2 {s^jyi [2 Crnk {$U^) nra 

+2 

' ru n 

(du) km + (dU) (dU) nm } 



15(q-1 + 2/x) 

a* (g — 1 + A 1 ) 3 
35 (q- 1 + 3//) (9- l) 2 



[l-/i + e(l + 2/i)] , 

(A.14) 
[2 + 3 (2 - 3/i) sin 2 0] , 
(A.15) 



is the angle between the unit vector e and the back- 
ground vorticity W due to the large-scale shear. Equa- 
tions (|A~T4|) - (|XT5)) yield Eq. (S) given in Sect. III. 
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